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Abstract: In the competing risks problem, an important role is played by the 
cumulative incidence function (GIF), whose value at time t is the probability 
of failure by time t from a particular type of failure in the presence of other 
risks. In some cases there are reasons to believe that the CIFs due to various 
types of failure are linearly ordered. El Barmi et al. [3] studied the estimation 
and inference procedures under this ordering when there are only two causes 
of failure. In this paper we extend the results to the case of k CIFs, where 
fc > 3. Although the analyses are more challenging, we show that most of the 
results in the 2-sample case carry over to this fc-sample case. 

1. Introduction 

In the competing risks model, a unit or subject is exposed to several risks at the 
same time, but the actual failure (or death) is attributed to exactly one cause. 
Suppose that there are fc > 3 risks and we observe (T, S) , where T is the time of 
failure and {6 = j} is the event that the failure was due to cause j, j = 1,2, . . . , k. 
Let F be the distribution function (DF) of T, assumed to be continuous, and let 
5 = 1 — F be its survival function (SF). 

The cumulative incidence function (GIF) due to cause j is a sub-distribution 
function (SDF), defined by 

(1.1) F,{t) = P[T<t, S = j], j = l,2,...,k, 

with F{t) — J2j The cause specific hazard rate due to cause j is defined by 

Xjit) = lim^ ^P[t <T<t + At, 6^ j\T>t], 1,2,..., k, 
and the overall hazard rate is X{t) = J2j ^ji^)- The GIF, Fj{t), may be written as 

(1.2) Fj{t) = [ Xj{u)S{u)du. 

Jo 

Experience and empirical evidence indicate that in some cases the cause specific 
hazard rates or the GIFs are ordered, i.e., 

Al < A2 < • • • < Afe or Fi < F2 < ■■■ < Fk. 
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The hazard rate ordering impHes the stochastic ordering of the CIFs, but not vice 
versa. Thus, the stochastic ordering of the CIFs is a milder assumption. El Barmi et 
al. [3] discussed the motivation for studying the restricted estimation using several 
real life examples and developed statistical inference procedures under this stochas- 
tic ordering, but only for k = 2. They also discussed the literature on this subject 
extensively. They found that there were substantial improvements by using the re- 
stricted estimators. In particular, the asymptotic mean squared error (AMSE) is 
reduced at points where two CIFs cross. For two stochastically ordered DFs with 
(small) independent samples, Rojo and Ma [17] showed essentially a uniform reduc- 
tion of MSE when an estimator similar to ours is used in place of the nonparametric 
maximum likelihood estimator (NPMLE) using simulations. Rojo and Ma [17] also 
proved that the estimator is better in risk for many loss functions than the NPMLE 
in the one-sample problem and a simulation study suggests that this result extends 
to the 2-sample case. The purpose of this paper is to extend the results of El Barmi 
et al. [3] to the case where fc > 3. The NPMLEs for k continuous DFs or SDFs under 
stochastic ordering are not known. Hogg [7] proposed a pointwise isotonic estimator 
that was used by El Barmi and Mukerjee [4] for k stochastically ordered continuous 
DFs. We use the same estimator for our problem. As far as we are aware, there are 
no other estimators in the literature for these problems. In Section 2 we describe 
our estimators and show that they are strongly uniformly consistent. In Section 3 
we study the weak convergence of the resulting processes. In Section 4 we show that 
confidence intervals using the restricted estimators instead of the empiricals could 
possibly increase the coverage probability. In Section 5 we compare asymptotic bias 
and mean squared error of the restricted estimators with those of the unrestricted 
ones, and develop procedures for computing confidence intervals. In Section 6 we 
provide a test for testing equality of the CIFs against the alternative that they are 
ordered. In Section 7 we extend our results to the censoring case. Here, the results 
essentially parallel those in the uncensored case using the Kaplan-Meier [9] estima- 
tors for the survival functions instead of the empiricals. In Section 8 we present an 
example to illustrate our results. We make some concluding remarks in Section 9. 

2. Estimators and consistency 

Suppose that we have n items exposed to k risks and we observe (Tj, Si), the time 
and cause of failure of the ith item, 1 < i < n. On the basis of this data, we wish 
to estimate the CIFs, Fi, F2, . . . , Fk, defined by (|1.1|) or (|1.2p . subject to the order 
restriction 



It is well known that the NPMLE in the unrestricted case when fc = 2 is given by 
(see Peterson, [12]) 



and this result extends easily to k > 2. Unfortunately, these estimators are not guar- 
anteed to satisfy the order constraint (|2.ip . Thus, it is desirable to have estimators 
that satisfy this order restriction. Our estimation procedure is as follows. 

For each t, define the vector F{t) = (A(t), hit), Fk{t))'^ and let I = {x e 
RI^ : Xi < X2 < ■ ■ ■ < Xk}, a closed, convex cone in i?*^. Let i?(x|X) denote the least 



(2.1) 



Fi<F2<---<Fk. 



(2.2) 
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squares projection of x onto T with equal weights, and let 



A?;[F; r, s] — 



s — r + 1 
Our restricted estimator of Fi is 

(2.3) F* =-ma^-{mnAv[f-r,s\=E{{Fi,...,Fkf\I),, 1 < i < k. 

r<i s>i 

Note that for each t, equation (|2.3p defines the isotonic regression of {Fi{t)}^^^ 
with respect to the simple order with equal weights. Robertson et al. [13] has a 
comprehensive treatment of the properties of isotonic regression. It can be easily 
verified that the F*s are CIFs for all i, and that X^Li ^ti^) = where F is the 
empirical distribution function of T, given by F{t) = ^i'^i — 0/*^ ^■ 

Corollary B, page 42, of Robertson et al. [13] implies that 

max \F*(t) - F,(t)\ < max \Fi(t) - F,(t)\ for each 
i<j<fc i<i<fc 

Therefore || F* — Fi \\< maxi<j<fe||/'j — Fj\\ for all i where ||.|| is used to denote 
the sup norm. Since — ^0 a.s. for all i, we have 

Theorem 2.1. P[||-F;* - as n ^ oo, i = 1,2, . . . , k] = 1. 

If fc = 2, the restricted estimators of Fi and F2 are Fj* = Fi A F/2 and F2 = 
Fi V F/2, respectively. Here A (V) is used to denote max (min). This case has been 
studied in detail in El Barmi et al. [3]. 



3. Weak convergence 

Weak convergence of the process resulting from an estimator similar to (|2.3p when 
estimating two stochastically ordered distributions with independent samples was 
studied by Rojo [15]. Rojo [16] also studied the same problem using the estimator 
in (|2.3p . Praestgaard and Huang [14] derived the weak convergence of the NPMLE. 
El Barmi et al. [3] studied the weak convergence of two CIFs using (|2.3p . Here we 
extend their results to the fc-sample case. Define 

Z.,n = MF^ - F,] and Z*„ = V^[F* ~ F,], i ^ 1,2, . . . , k. 

It is well known that 

(3.1) {Z In, Z2n, ■ ■ ■ , Zkn) =^ {Zl, Z2, ■ ■ ■ , Zk) , 

a fc-variate Gaussian process with the covariance function given by 

CoviZ,{s), Zj{t)) = F,{s)[S,j - Fj{t)], l<i,j<k, for s < t, 

where Sij is the Kronecker delta. Therefore, Zi = Bf{Fi) for all i, the B^s being 
dependent standard Brownian bridges. 

Weak convergence of the starred processes is a direct consequence of this and the 
continuous mapping theorem. First, we consider the convergence in distribution at 
a fixed point, t. Let 



(3.2) 



Su^{j:Fj{t) = F,{t)}, z = l,2,...,fc. 
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Note that Su is an interval of consecutive integers from {1, 2, . . . , k}, Fj{t) — Fi{t) = 
for j e Sit, and, as n ^ oo, 

(3.3) Vn[Fj{t) - Fi{t)] oo, and V7i[Fj{t)- F,{t)] -oo, 

for j > i*{t) and j < respectively, where = min{j ; j G Sn} and 

i*{t) = max{j -.jeSit}. 

Theorem 3.1. Assume that (j2.ip holds and t is fixed. Then 

(ZLit), z^t), z*,,Mf {zm,zi{t), z*,{t)f, 

where 

/n A\ ry*l,\ • ^ {r<j<s} Z jit) 

(3.4) z, (t) = max mm — ^ ~ -* . 

i.(t)<r<ii<s<i*(t) s-r+l 

Except for the order restriction, there are no restrictions on the Fi'& for the conver- 
gence in distribution at a point in Theorem 2. For A; = 2, if the Fi% are distribution 
functions and the FiS are the empiricals based on independent random samples of 
sizes ni and 712, then, using restricted estimators F*s that are slightly different from 
those in (|2.3p . Rojo [15] showed that the weak convergence of (ZJ'„^,Z|„^) fails if 
Fi(6) = F2(6) and Fi < F2 on {b, c] for some b<c with < ^2(6) < ^2(0) < 1. El 
Barmi et al. [3] showed that the same is true for two CIFs. They also showed that, 
if Fi < F2 on (0, b) and Fi = F2 on [b, 00), with Fi{b) > 0, then weak convergence 
holds, but the limiting process is discontinuous at b with positive probability. Thus, 
some restrictions are needed for weak convergence of the starred processes. 

Let a {di) be the left (right) endpoint of the support of Fi, and let Si = {j : Fj = 
Fi} for i = 1, 2, . . . , A:. In most applications Ci = 0. Letting i* = max{j : j € Si}, 
we assume that, for i = 1, 2, . . . , fc — 1, 

(3.5) inf [Fj(i) -Fi(i)] > for all > and j > r. 

Ci+r]<t<di—ri 

Note that i E Si for all i. Assumption (|3.5p guarantees that, if Fj > Fi, then, either 
Fj = Fi or Fj{t) > Fi{t), except possibly at the endpoints of their supports. This 
guarantees that the pathology of nonconvergence described in Rojo [15] does not 
occur. Also, from the results in El Barmi et al. [3] discussed above, if di = dj for 
some i ^ j ^ Si, then weak convergence will hold, but the paths will have jumps at 
di with positive probability. We now state these results in the following theorem. 

Theorem 3.2. Assume that condition (j2.ip and assumption (13. 5p hold. Then 



where 



l^lriJ '^2riJ ■ • • ' ^fcni =^ K'^ 1 y -^2 ^ ■ ■ ■ ^ k) y 
-7* • ^{r<j<s} Zj 

Z: = max mm ^ j 



,<r<ii<s<i* S — r + 1 

Note that, if Si = {i}, then Z*„ ^> Zi under the conditions of the theorem. 
4. A stochastic dominance result 

In the 2-sample case. El Barmi et al. [3] showed that \Z* \ is stochastically dominated 
by \Zj\ in the sense that 

P[\Z*{t)\ <u]> P[\Zj{t)\ <u], j = 1,2, for aU u > and for all t. 
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if < Fi{t) = F2{t) < 1. This is an extension of Kelly's [10] result for indepen- 
dent samples case, but restricted to fc = 2; Kelly called this result a reduction of 
stochastic loss by isotonization. Kelly's [10] proof was inductive. For the 2-sample 
case, El Barmi et al. [3] gave a constructive proof that showed the fact that the 
stochastic dominance result given above holds even when the order restriction is 
violated along some contiguous alternatives. We have been unable to provide such 
a constructive proof for the fc-sample case; however, we have been able to extend 
Kelly's [10] result to our (special) dependent case. 

Theorem 4.1. Suppose that for some 1 < i < k, Su, as defined in (|3.2p . contains 
more than one element for some t with < Fi{t) < 1. Then, under the conditions 
of Theorem 3, 

P[\Z*{t)\ <u]> P[\Z,{t)\ < u] for all u > 0. 

Without loss of generality, assume that Su = {j : Fj{t) = Fi{t)} = {1, 2, . . . ,1} 
for some 2 < I < k. Note that {Zi{t)} is a multivariate normal with mean 0, and 

(4.1) Cov{Z,{t),Z,{t)) = Fi{t)[6,, - F^{t)l 1 < I,] < I. 

Also note that {Z*{t); 1 < j < fc} is the isotonic regression of {Zj{t) : 1 < j < fc} 
with equal weights from its form in (j3.4p . Define 

(4.2) X(')(t) = (Zi(t) - Z,{t), Z2{i) - Z,{i), . . . , Zi{t) - Z,it)f. 
Kelly [10] shows that, on the set {Z,{t) ^ Z*(t)}, 

(4.3) P[|Z*(0| < u|X«] > P[|Z,(t)| < u|X«] a.s. Vu > 0, 

using the key result that X'*^ {t) and Av{Z{t); 1, fc) are independent when the Zi{tys 
are independent. Although the Zi(t)s are not independent in our case, they are 
exchangeable random variables from (|4.ip . Computing the covariances, it easy to 
see that X'^'^(t) and Av{Z{t);l,k) are independent in our case also. The rest of 
Kelly's [10] proof consists of showing that the left hand side of (|4.3[) is of the form 
^(a + v)~ $(a — v), while the right hand side of ()4.3p is ^{b + v) — $(6 — v) using 
()4.2|) . where <i> is the standard normal DF, and b is further away from than a. 
This part of the argument depends only on properties of isotonic regression, and it 
is identical in our case. This concludes the proof of the theorem. 

5. Asymptotic bias, MSE, confidence intervals 

If Sit — {*} for some i and t, then Z*{t) = Zi{t) from Theorem 2, and they have 
the same asymptotic bias and AMSE. If Su has more than one element, then, 
for fc = 2, El Barmi et al. [3] computed the exact asymptotic bias and AMSE 
of Z*{t), i = 1,2, using the representations, Z^ = Zi + A {Z2 - Zi)/2 and 
Z2 = ^2 - A (Z2 - Zi)/2. The form of Z* in (13. 4p makes these computations 
intractable. However, from Theorem 4, we can conclude that E[Z*{t)]'^ < E[Zi{t)]'^, 
implying an improvement in AMSE when the restricted estimators are used. 

From Theorem 4.1 it is clear that confidence intervals using the restricted es- 
timators will be more conservative than those using the empiricals. Although we 
believe that the same will be true for confidence bands, we have not been able to 
prove it. 
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The confidence bands could always be improved by the following consideration. 
The 100(1 — a)% simultaneous confidence bands, [Li, Ui], for Fi, 1 < i < k, in the 
unrestricted case obey the following probability inequality 

P{F,e[L„U^]:l<i<k)>l~a. 

Under our model, Fi < F2 < ■ ■ ■ < Fk, this probability is not reduced if we replace 
[L,,U,] by [L*,f/*], where 

L* = max{Lj : 1 < j < i} and U* — min{C/j '■ i < j ^ k}, 1 < i < k. 
6. Hypotheses testing 

Let Hq : Fi — F2 ^ ■ ■ ■ = Fk and iJa : Fi < F2 <•••< . In this section we 
propose an asymptotic test of Hq against Ha — Hq. This problem has already been 
considered by El Barmi et al. [3] when k = 2, and the test statistic they proposed 
is T„ = ^/nsnp^yQ[F2{x) — Fi{x)]. They showed that under Hq, 

(6.1) lim P{T,a >t) = 2(1 - <I>(t)), t > 0, 

n — >oo 

where $ is the standard normal distribution function. 

For k > 2, we use an extension of the sequential testing procedure in Hogg [7] 
for testing equality of distribution functions based on independent random samples. 
For testing Hoj : Fi = F2 = ■ ■ ■ = Fj against Haj — Hoj , where Haj : Fi = F2 = 
■ ■ ■ — Fj^i < Fj,] — 2,3, . . . , k, we use the test statistic sup^>o Tjn{x) where 

T,n = V^V^[F, - Av[F; l,j - 1]], 

with Cj — k{j — We reject H^j for large values Tj„, that may be also written 
as 

T-jn = ^[Zjn - Av{Zn; l,j ~ 1)], 

where Z„ — {Z\n, Z2m ■ ■ ■ , ZknY" ■ By the weak convergence result in p.ip and the 
continuous mapping theorem, {T2n, T3„, . . . , Tkn)^ converges weakly to (T2, Ta, . . . , 
Tfc)^, where 

T,=^,[Z,-Av[Z;l,j-l]]. 

A calculation of the covariances shows that the T^'s are independent. Also note 
that 

T,^B,{F), 2<3<k, 

where the Bj's are independent standard Brownian motions and F = 2_/i=i = 
kFi under Hq. We define our test statistic for the overall test of Hq against Ha — Hq 

by 

T„ — max sup Tjn{x). 

2<J<*: x>0 

By the continuous mapping theorem, r„ converges in distribution to T, where 

T= max supTj(a;). 

2<i<fc x>0 
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Using the distribution of the maximum of a Brownian motion on [0, 1] (BiUingsley 
[2]), and using the independence of the -B^'s, the distribution of T is given by 

P{T >t) = l- P{snpTj{x) <t,j = 2,...,k) 

X 

k 

= l-\{P{snpB,{F{x))<t) 

= l-[2^{t)-lf-^. 
This allows us to compute the p- value for an asymptotic test. 

7. Censored case 

The case when there is censoring in addition to the competing risks is considered 
next. It is important that the censoring mechanism, that may be a combination 
of other competing risks, be independent of the k risks of interest; otherwise, the 
CIFs cannot be estimated nonparametrically. We now denote the causes of failure 
as (5 = 0, 1, 2, . . . , fc, where {5 — 0} is the event that the observation was censored. 

Let Ci denote the censoring time, assumed continuous, for the ith subject, and 
let Li —Tif\Ci. We assume that C^s are identically and independently distributed 
(IID) with survival function, 5c, and are independent of the life distributions, {Ti}. 
For the ith subject we observe (L^, Si), the time and cause of the failure. Here the 
{Li} are IID by assumption. 



7.1. The estimators and consistency 

For i — 1, 2, . . . , fc, let be the cumulative hazard function for risk j, and let 
A = Ai + A2 + • • • + Afc be the cumulative hazard function of the life time T. For the 
censored case, the unrestricted estimators of the CIFs are the sample equivalents 
of (II. 2p using the Kaplan-Meier [9] estimator, S*, of 5 = 1 — i^: 

(7.1) F,{t)= [ S{u)dA,{u), j^l,2,...,k, 

Jo 

with F = Fi + F2 + ■ ■ ■ + Fk, where S is chosen to be the left-continuous version 
for technical reasons, and A^ is the Nelson-Aalen estimator (see, e.g., Fleming and 
Harrington, [5]) of Aj. Although our estimators use the Kaplan-Meier estimator of 
S rather than the empirical, we continue to use the same notation for the various 
estimators and related entities as in the uncensored case for notational simplicity. 
As in the uncensored case, we define our restricted estimator of Fi by 

max min Av [F ; r , s] 

r<i s>i 

EiiFi,...,Fkf\I)„ l<i<k. 
Let 

At) = P[L, >t]= P[T, > t, C, >t] = Sit)Scit). 

Strong uniform consistency of the i^*s on [0, b] for all b with 7r(6) > follows from 
those of the F^'s [see, e.g., Shorack and Wellner [18], page 306, and the corrections 



F* = 



(7.2) 
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posted on the website given in the reference] using the same arguments as in the 
proof of Theorem 2 in the uncensored case. 



7.2. Weak convergence 

Let Zjn = V^[F] - Fj] and Z*„ = \/n[F* - Fj], j = 1, 2, . . . , fc, be defined as in 
the uncensored case, except that the unresticted estimators have been obtained via 
(|7.ip . Fix b such that 7r(&) > 0. Using a counting process-martingale formulation, 
Lin [11] derived the following representation of Zj„ on [0,6]: 

t 



Y{u) '^Vo Y{u) 



F,iu)j:'dM,iu) 







Y{u) 



where 



Y{t) = ^ I{Lj > t) and M,{t) = ^ I{Lj < t) - ^ / I{Lj > u)dA,{u), 
j=i j=l 3=1 ''° 

the Mi's being independent martingales. Using this representation, El Barmi et al. 
[3] proved the weak convergence of (Zj„,Z2n)^ to a mean-zero Gaussian process, 
(Zi, Z2), with the covariances given in that paper. A generalization of their results 
yields the following theorem. 

Theorem 7.1. The process (Zi„, Z2„, . . . , ZknV ^ i^i, ^2, • • ■ , ^fc)^ on [0, b]'', 
where {Zi, Z2, ■ ■ ■ , Zk)^ is a mean-zero Gaussian process with the covariance func- 
tions, for s < t, 

(u))] '^^^(^ 

^ [\F.,{u)-F,{s)m{u)-F,{t)]- 



Cov{ZM,Zm = I [l-F,is)^Y.F^{umi^F,{t)-J2F,iu))]- 

- in) ■ 
and, for i^ j, 

(u) 



(7.4) 



Cov{Z,{s),Z,{t))^ f [l-F,{s)-J2Fi{u)][F,{u)-F,{t)] — 

^ \\f,{s)-F, 



/ yr,ys^-r,{u)\{F,{t)~F^(u)\'^^'^'''^ 



7r(u) 



The proofs of the weak convergence results for the starred processes in Theorems 
3.1 and 3.2 use only the weak convergence of the unrestricted processes and iso- 
tonization of the estimators; in particular, they do not depend on the distribution 
of {Z\, . . . , Zk)^ . Thus, the proof of the following theorem is essentially identical to 
that used in proving Theorems 3.1 and 3.2; the only diS'erence is that the domain 
has been restricted to [0,6]'^. 
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Theorem 7.2. The conclusions of Theorems 3.1 and 3.2 hold for (2'J'^,Z2„, 
^kn)'^ defined above on [0,6]'^ under the assumptions of these theorems. 

7.3. Asymptotic properties 

In the uncensored case, for a f > and an i such that < < 1, if Su = 

{!,...,/} and l>2, then it was shown in Theorem 4 that 

P{\Z*{t)\ <u)> P{\Zi{t)\ < u) for all u>0. 

The proof only required that {Zj{t)} be a multivariate normal and that the ran- 
dom variables, {Zj{t) : j £ Su}, be exchangeable, which imply the independence of 
XW(t) and AviZ{t); 1,1), as defined there. Noting that Fj{t) = Fi{t) for all j € S^t, 
the covariance formulas given in Theorc^ni 7.1 show that the multivariate normal- 
ity and the exchangeability conditions hold for the censored case also. Thus, the 
conclusions of Theorem 4.1 continue to hold in the censored case. 

All comments and conclusions about asymptotic bias and AMSE in the uncen- 
sored case continue to hold in the censored case in view of the results above. 



7.4- Hypothesis test 

Consider testing Hq : Fx = = ■ ■ ■ = against Ha — Hq, where Ha ■ Fi < F2 < 
■ ■ ■ < Fk, using censored observations. As in the uncensored case, it is natural to 
reject Hq for large values of T„ = max2<j<fe sup^^QTjn{x), where 

Tjn{x) = V^^[Fj{x) - Av{F{xy, l,j - 1)] 
= ^[Zjn{x) - Av{Zn{xy, 1,3 - 1)] 

with Cj = k{j — is used to test the sub-hypothesis Hqj against Haj — Hqj, 2 < 
j < k, asm the uncensored case. Using a similar argument as in the uncensored case, 
under Hq, (T2„, T^n, ■ ■ ■ , Tkn) converges weakly (T2, T3, . . . , Tfc)^ on [0, 6]^=, where 
the Tj's are independent mean zero Gaussian processes. For s < t, Cov{Ti{s) ,Ti{t)) 
simplifies to exactly the same form as in the 2-sample case in El Barmi et al. [3] : 



Cov{Ti{s),Ti{t)) = I S{u 
Jo 

The limiting distribution of 



dA{u) 



Tn = max sup[Zjn{x) - Av{Zn{x); l,j - 1)] 

2<j<fcx>0 

is intractable. As in the 2-sample case, we utilize the strong uniform convergence 
of the Kaplan-Meier estimator, Sc, of Sc, to define 

T*^{t) = V^^ I JSc{u) d[Fj{x) - Av{¥{x), l,j - = 2, 3, . . . , fc, 
Jo 

and define T* = max2<j<kSup^^QT*^{x) to be the test statistic for testing the 
overall hypothesis of Ho against Ha — Hq- By arguments similar to those used in 

the uncensored case, {T^^, r{„, . . . , T^J"^ covcrgcs weakly to {T^, T3*, . . . , T^)"^, a 
mean zero Gaussian process with independent components with 

T*^Bj{F), 2<j<k, 
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where Bj is a standard Brownian motion, and T* converges in distribution to a 
random variable T*. Since Tj here and Tj in the uncensored case (Section 6) have 
the same distribution, 2 < j < k, T* has the same distribution as T in Section 6, 
i.e., 

P{T* >t) = l- [2$(i) - 1]*=-^ 

Thus the testing problem is identical to that in the uncensored case, with r„ of 
Section 6 changed to T* as defined above. This is the same test developed by Aly 
et al. [1], but using a different approach. 



8. Example 



We analyze a set of mortality data provided by Dr. H. E. Walburg, Jr. of the 
Oak Ridge National Laboratory and reported by Hoel [6]. The data wcire obtained 
from a laboratory experiment on 82 RFM strain male mice who had received a 
radiation dose of 300 rads at 5-6 weeks of age, and were kept in a conventional 
laboratory environment. After autopsy, the causes of death were classified as thymic 
lymphoma, reticulum cell sarcoma, and other causes. Since mice are known to be 
highly susceptible to sarcoma when irradiated (Kamisaku et al [8]), we illustrate our 
procedure for the uncensored case considering "other causes" as cause 2, reticulum 
cell sarcoma as cause 3, and thymic lymphoma as cause 1, making the assumption 
that -Fi < F2 < F3. The unrestricted estimators are displayed in Figure 1, the 
restricted estimators are displayed in Figure 2. We also considered the large sample 
test of Hq : Fi = F2 = F3 against Ha — Hq, where Ha : Fi < F2 < F^, using the 
test described in Section 6. The value of the test statistic is 3.592 corresponding to 
a p-value of 0.00066. 



9. Conclusion 



In this paper we have provided estimators of the GIFs of k competing risks under 
a stochasting ordering constraint, with and without censoring, thus extending the 




Fig 1. Unrestricted estimators of the cumulative incidence functions. 
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Fig 2. Restricted estimators of the cumulative incidence functions. 



results for fc = 2 in El Barmi et al. [3]. We have shown that the estimators are 
uniformly strongly consistent. The weak convergence of the estimators has been 
derived. We have shown that asymptotic confidence intervals are more conservative 
when the restricted estimators are used in place of the empiricals. We conjecture 
that the same is true for asymptotic confidence bands, although we have not been 
able to prove it. We have provided asymptotic tests for equality of the CIFs against 
the ordered alternative. The estimators and the test are illustrated using a set of 
mortality data reported by Hoel [6] . 
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